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pg ' Abstract 

^ ' Preliminary results of a systematic study on the simultaneous optical-to- X-rays 

variability in blazars are presented. Data from Swift observations of four bright 
7-ray blazars (3C 279, ON +231, S5 0716+71, PKS 2155-304) have been analyzed, 
compared, and discussed. Specifically, 3C 279 shows a variable flattening in the low 
1-^ ■ energy part of the X-ray spectrum that appears to be confined in a specific X-ray 

vs optical/UV fluxes region. Some implications are shortly analyzed. 
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S ■ Active galactic nuclei with relativistic jets pointed to or within a few degrees 

O ■ of the observer's direction are generally named blazars. However, this name 

collects a variety of phenomena, most of them are still to be understood. The 
k> ; spectral energy distribution (SED - u vs uFy) of these sources always dis- 

^ ■ plays a specific shape, which is the result of two types of emission: the first 



one, extending from radio to soft X-rays, is thought to be due to the syn- 
chrotron emission from relativistic electrons; while the second one, peaking 
at energies in the X-/7— ray band, is generally attributed to inverse- Compton 
(IC) emission, resulting from the transfer of energy from the same popula- 
tion of relativistic electrons impinging on low-energy seed photons. The latter 
population of photons can be the same synchrotron radiation (Synchrotron 
Self-Compton, SSC) or a source external to the jet, like the accretion disk 
or the broad-line region (External Compton, EC). Fossati et al. (1998) and 
Ghisellini et al. (1998) suggested that blazars follow a trend - the hlazar se- 
quence - where sources with low bolometric luminosity have both peaks at 
high frequencies (soft X-rays for synchrotron and TeV for IC), while, as the 
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luminosity increases, both peaks shift to lower frequencies (optical/infrared 
for synchrotron and MeV/GeV for IC). See, however, Padovani (2007) and 
references therein for a discussion. Recent updates of the blazar sequence can 
be found in Ghisellini & Tavecchio (2008) and Maraschi et al. (2008). 

The shape and normalization of the SED during quiescence can be differ- 
ent when the source is in outburst and it is still to be understood how and 
when these changes occur. To investigate this problem, I started a study of 
the simultaneous optical/UV/X-rays variability of a sample of 7— ray blazars 
(detected both at MeV/GeV and TeV energies) to search for some common 
patterns of variability. Here I present the preliminary results. 



2 Sample selection and analysis 



Blazars are generally studied with multiwavelength (MW) campaigns, which 
are triggered by events of outburst and therefore the quiescence state is often 
neglected. However, to understand blazars, it is necessary to know how the 
source behaves during all its lifetime. Therefore, I focus on bright blazars, 
which can be almost always observed by Swift satellite (Gehrels et al. 2004) 
both during quiescence and outburst. It is also necessary to select blazars 
observed many times (to have a sufficient number of points) and with the 
same instruments settings (to have the same type of MW information). For 
the moment, I found sufficient data for the 4 blazars here shown, but work is 
in progress to enlarge the sample. 

Swift data have been processed, cleaned and analyzed by using the software 
package HEASof t 6 . 5 and the calibration data base CALDB updated on June 
25, 2008. Data from the X-Ray Telescope (XRT, 0.2 - 10 keV, Burrows et al. 
2005) were analyzed with the xrtpipeline task with standard parameters. 
XRT was set in photon counting (PC) mode for the observations of 3C 279, 
ON +231, S5 0716 + 714 and single to quadruple pixel events were selected 
(grades — 12). In the case of PKS 2155 — 304, XRT was set in window tim- 
ing (WT) mode for most of the pointings and events with grades — 2 were 
considered. A small numbers of observations in PC mode were performed: in 
this case, being the source count rate so high to be affected by pile-up, it was 
necessary to apply a correction by removing the central region of the point- 
spread function. Therefore, the source extraction region was an annulus with 
10" and 50" of internal and external radius, respectively. The extracted spec- 
tra were rebinned to have at least 20 — 30 counts per energy bin, depending 
on the available statistics, and fitted with models in xspec v. 11.3.2ag. I 
used a single redshifted power-law or a broken power-law models, with Galac- 
tic absorption from Kalberla et al. (2005). Sometimes, the low energy part of 
the X-ray spectrum (generally < 2 keV) shows a curvature due to a deficit of 



Table 1 

Summary of Swift /XKT observations of 3C 279 (z = 0.5362, Afn.Gal = 2.05 X 10^0 cm^^^ Kalberla ct al. 
2005). Columns: (1) Date of observation [DD-MM-YYYY]; (2) Photon index (for the power-law model) 
or soft photon index (for the broken power-law model); (3) Break energy £"(, for the broken power-law 
model or folding energy Ef of the exponential roll-off for the single power-law model plus low energy 
photon deficit [keV]; (4) Hard photon index for the broken power-law model; (5) Model normalization 
[10^'* ph cm^'^ s^-"^ kcV^^]; (6) Reduced x^ and degrees of freedom of the spectral fit; (7) Notes. 
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photons that cannot be fitted with a broken power-law model. In this case, 
the single power-law model is multiplied by an exponential roll-off in the form 
exp{—Ej/E'^), where Ef is the folding energy of the roll-off. The fits are sum- 
marized in the Tables 1, 2, 3 and 4. 



-5 -"5 



Optical and ultraviolet data of UVOT (Roming et al. 2005) in the filters 
B (4329 A) and UVWl (2634 A) were considered. The snapshots of every 
individual pointing were integrated with uvotimsum task and then analyzed 
with the uvotsource task. The source region was a circle with radius of 5" 
for B filter and 10" for UVWl filter. The background was extracted from an 
annular region centered on the blazar, with internal radius of 7" for optical 
filters and 12" for UV. The external radius was selected on the basis of the 
presence or not of nearby contaminating sources and was 40" for ON -1-231, 
50" for 3C 279 and PKS 2155 - 304, 60" for S5 0716 + 714. 

To emphasize the source variability and minimize the systematic or contami- 



Table 2 

Summary of Swift/XKT observations of ON +231 {z = 0.102, Afn.Gal = 2.04 X 10^0 cm^^^ Kalbcrla et 
al. 2005). Columns: (1) Date of observation [DD-MM-YYYY]; (2) Photon index (for the power-law model) 
or soft photon index (for the broken power-law model); (3) Break energy Ei, for the broken power-law 
model or folding energy Ef of the exponential roll-off for the single power-law model plus low energy 
photon deficit [keV]; (4) Hard photon index for the broken power-law model; (5) Model normalization 
[lO^'* ph cm^'^ s^i kcV^i]; (6) Reduced x^ and degrees of freedom of the spectral fit; (7) Notes. 
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Table 3 

Summary of Swift/XKT observations of 85 0716 + 714 (z = 0.31, A^H.Gal = 3.11 X 10^° cm-^, Kalberla et 
al. 2005). Columns: (1) Date of observation [DD-MM-YYYY]; (2) Photon index (for the power-law model) 
or soft photon index (for the broken power-law model); (3) Break energy £"(, for the broken power-law 
model or folding energy Ef of the exponential roll-off for the single power-law model plus low energy 
photon deficit [keV]; (4) Hard photon index for the broken power-law model; (5) Model normalization 
[10~^ ph cm~^ s~l keV~l]; (6) Reduced x^ and degrees of freedom of the spectral fit; (7) Notes. 
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nating effects, the selected parameters (X-ray, UV and Optical flux densities. 



Table 4 

Summary of Swift/XRT observations of PKS 2155 - 304 {z = 0.116, iVH.Gal = 1-48 X lO^" cm-2, Kalbcrla 
et al. 2005). Columns: (1) Date of observation [DD-MM-YYYY]; (2) Photon index (for the power-law 
model) or soft photon index (for the broken power-law model); (3) Break energy _E(, for the broken power- 
law model or folding energy E-f of the exponential roll-off for the single power-law model plus low energy 
photon deficit [keV]; (4) Hard photon index for the broken power-law model; (5) Model normalization 
[10^^ ph cm^^ s^^ keV^-'^]; (6) Reduced x^ and degrees of freedom of the spectral fit; (7) Notes. 
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r) have been normalized to their average value (indicated between <>). The 
results are shown in Fig. 1, 2, 3 and 4, where red points refer to X-ray spectra 
fitted with a single power-law model and black points indicate data fitted with 
a broken power-law model or with a single power-law with low-energy expo- 
nential roll-off. In the case of a broken power-law model, the photon index 
used in the figure is the soft one, i.e. below the energy break. 
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Fig. 1. X-ray flux density at 1 keV versus T (top panel), UV flux density at 2634 A{center panel), and 
optical flux density at 4329 A{bottom panel) for 3C 279. 
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Fig. 2. X-ray flux density at 1 keV versus T (top panel), UV flux density at 2634 A{center panel), and 
optical flux density at 4329 A{bottom panel) for ON +231. 
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Fig. 3. X-ray flux density at 1 keV versus T (top panel), UV flux density at 2634 A{center panel), and 
optical flux density at 4329 A{bottom panel) for S5 0716 + 71. 



> 



<r>=2.6±0.1 

<Fx>=(3.9±0.2)xlO-- ph cm-= s"' keV"' 




iBiHSif ^ 



Giant TeV flare (200S) 



0.5 1 1.5 

Normalized Flux Density at 1 keV 



<Fm,>=17.4±0.4mJy 
<F^>=(3.9±0.2)xlO-2 ph cm-' s"' keV"' 



1 |i , — h- 



+ 



0.5 1 1.5 

Normalized Flux Density at I keV 



<Fu>=25.43±0.09 mjy 
<Fx>=(3.9±0.2)xlO-- ph cm-= s"' keV"' 



0.5 I 1.5 

Normalized Flux Density at I keV 



Fig. 4. X-ray flux density at 1 keV versus T (top panel), UV flux density at 2634 A{center panel), and 
optical flux density at 4329 k(bottom panel) for PKS 2155 - 304. 



3 Discussion 



As expected, all the sources display correlated variability from optical to X- 
ray frequencies, since the radiation emitted by a blazar is dominated by the 
relativistic jet component, but there are some specific points worth noting. 

The most interesting case is 3C 279, whose observations were performed during 
two MW campaigns, one in mid- January 2006 (Bottcher et al. 2007) triggered 
by high optical state {R = 14.5), while the other was performed in January 
2007, with an even higher optical state {R ~ 13). The black points in Fig. 1, 
left column, indicate the deficit of photons in the low energy part of the X- 
ray spectrum, which appears to be confined in a "branch" separated by the 
remaining. This photon deficit can be modeled either with a broken power- 
law with Tsoft < ^hard (low-euergy spectral flattening) and a break between 
~ 1.5 — 1.9 keV, or with a single power-law with low-energy exponential roll-off 
and folding energy of ~ 0.35 keV. These changes occur on timescale of days. 
A similar effect, i.e. changes in the low-energy spectral flattening, has been 
already noted by Bianchin et al. (2008) in the case of the blazar PKS 2149—306 
{z = 2.345). 

The nature of this flattening, often observed in high-redshift blazars, has been 
examined and is still under investigation. The low-energy spectrum {E < 
2 keV) can give important information about the very low end of the electron 
distribution, since - according to the EC models - a break is expected at 
t'br,obs - ^'BBrbuik7min/(l + z) , where Tbuik is the bulk Lorentz factor, 7min 
is the minimum Lorentz factor of the electron distribution, ubb is the peak 
frequency of the emission of the accretion disk (modeled as a black body). In 
the case of 3C 279 and many other fiat-spectrum radio quasars (FSRQ), the 
typical values of the parameters are (e.g. Ballo et al. 2002): Fbuik = 10 — 20, 
7mm = 1) ^BB ~ 10^^ Hz, z = 0.5362, thus resulting in z/br,obs ~ 10^^ Hz. A 
more detailed discussion of the low-energy photon deficit in the framework of 
a natural curvature of the blazar continuum can be found in Tavecchio et al. 
(2007) and Ghisellini et al. (2007). 

In reality, things are not so easy and self-explaining: indeed, at these energies 
there could be contributions from SSC (as in the case of 3C 279, Ballo et al. 
2002) or from other processes, like the bulk motion Comptonization (Celotti 
et al. 2007). In addition, since this break appears in the X-ray spectrum as a 
photon deficit at low-energies, some authors suggested that it can be explained 
as an absorption (cold or ionized) in addition to the Galactic column and 
intrinsic to the source (e.g. Cappi et al. 1997, Fiore et al. 1998, Piconcelli & 
Guainazzi 2005, Page et al. 2005, Yuan et al. 2006). 

The data on 3C 279 shown in Fig. 1, left column, favor the hypothesis of 
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the natural break in the SED, which is the curvature of the continuum of a 
typical FSRQ. The points of the flattening branch (Fig. 1) are correlated each 
other at all the analyzed frequencies (optical to X-rays), as expected from the 
emission of the relativistic jet. Would it be something external to the jet, like 
an absorption intrinsic to the source, it is expected to find no correlation with 
the changes in the continuum of the source. 

Also the black point in ON +231 (Fig. 1, right column) indicates a flattening at 
low energies in the X-ray spectrum and refers to an observation performed on 
March 14, 2008, just one day after the detection of a VHE flare by VERITAS 
(Swordy et al. 2008). However, since ON +231 is an intermediate blazar, it 
is not expected to have the same behavior of 3C 279. One single point is not 
sufficient for a fruitful investigation, but surely the fact that the low-energy 
flattening occurred in proximity of a VHE flare, make it an interesting case 
worth studying. 

On the other hand, the black points in S5 0716 + 71 (Fig. 4, left column) 
indicate flt with a broken power-law model, with Tsoft > ^hard, where the 
soft component is due to the tail of the synchrotron emission and the hard 
X-ray emission is due to the emergence of the inverse- Compt on component 
(e.g. Foschini et al. 2006). The points at highest X-ray flux were measured a 
few days after the VHE detection by MAGIC (Teshima et al. 2008), while the 
observations during the AGILE pointings had lower X-ray fluxes (cf. Giommi 
et al. 2008). 

Interestingly, despite the fact that both ON +231 and S5 0716 + 71 are in- 
termediate blazars, the former was detected at TeV energies during a low 
optical/UV state, while the contrary occurred for the latter. 

PKS 2155 — 304 was extensively observed during the giant TeV flares in July 
2006 and a more detailed analysis has been presented elsewhere (Foschini et 
al. 2007). Here I conflrm the previous analyses. 



4 Conclusion 



I started a systematic study of the simultaneous optical-to- X-ray variability of 
blazars and presented here the preliminary results. The most interesting result 
concerns the FSRQ 3C 239, which shows a clear pattern when a low-energy 
photon deflcit appears in the X-ray spectrum. Although this deflcit has been 
interpreted in different ways, the data available here favor the hypothesis that 
is due to the natural curvature of the continuum in a typical FSRQ. 

Future steps of this research are to enlarge the sample and continue the mon- 
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itoring, taking advantage of the unique capabilities of the Swift satelhte. In 
addition, the successful beginning of the operations of the Fermi Gamma- 
ray Space Telescope during the summer of 2008, will give, in the near future, 
the possibility to include also the data at 7— rays, which are of paramount 
importance in defining the true state of activity of the source. 
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